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In this paper, we present analytical and numerical calculations of the 

7r resonance in the t-J model. We show in detail how the it resonance in 

the particle-particle channel couples to and appears in the dynamical spin 

correlation function in a superconducting state. The contribution of the ir 

resonance to the spin excitation spectrum can be estimated from general 

model-independent sum rules, and it agrees with our detailed calculations. 

The results are in overall agreement with the exact diagonalization studies of 

the t-J model. Earlier calculations predicted the correct doping dependence of 

the neutron resonance peak in the YBCO superconductor, and in this paper 

detailed energy and momentum dependence of the spin correlation function 

is presented. The microscopic equations of motion obtained within current 

formalism agree with that of the SO (5) nonlinear sigma model, where the 

7r resonance is interpreted as a pseudo Goldstone mode of the spontaneous 

SO (5) symmetry breaking. 

I. INTRODUCTION 

Of many fascinating experiments on high T c superconductors, the resonant neutron scat- 
tering peak observed in the YBCO family is an extremely striking one [3]-||. It was first 
observed in the optimally doped YBCO materials. The mode exists only in a narrow region 
in reciprocal space near (n/a, ir/b, tt /c), where a and b are the lattice constants in the Cu02 
plane and c is the distance between two neighboring Cu02 planes in a unit cell. (In the fol- 
lowing, we will set these lattice constants to unity to simplify notations). The energy of the 
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resonance is AlmeV and it disperses weakly in reciprocal space. Perhaps the most striking 
property of this mode is its disappearance above T c . More recently, this type of collective 
mode has also been observed in the underdoped families of the YBCO superconductors. 
Here the energy of this mode is 33meV and 25meV, for materials with T c values of Q2K 
and 52K respectively. While the mode energy decreases monotonically with T c , the mode 
intensity increases as T c decreases. Compared with the AlmeV peak, these modes also have 
a broader spectral distribution below T c . In these underdoped materials the resonance is 
also observed above T c where it becomes significantly broader. All the modes have been 
observed in the neutron spin flip channel, and more recently, the 41meV mode was seen to 
broaden under a uniform magnetic field ||, both indicating that the modes are spin triplets. 

These striking resonances have generated wide theoretical interests and a number of 
theoretical ideas have been suggested in order to explain their properties |j7Hl6fl. We believe 
that one key ingredient is the coupling of the neutron to the particle-particle (p-p) channel 
which occurs in the superconducting (SC) state via the condensate. In particular for a d x 2__ y 2 
gap, the coherence factor [1 — A k+q A k / E k E k+q }/2 for Q = (jr, tc, 7r) goes to unity at threshold 
rather than vanishing as it would for an s-wave superconductor P|,|TTfl. Furthermore, two of us 
argued that the p-p interactions in this channel leads to a sharp resonance which was called 
the 7r-mode [0. In the normal state, the resonance is decoupled from the neutron scattering, 
but can in principle be observed in pair tunneling experiments |Tj|. This theory predicted 



the doping dependence of the mode energy and intensity which was subsequently verified 
experimentally fl9f. This picture was also later verified in detailed numerical calculations 



of the Hubbard and the t-J models pO| , pT 



In this paper, we study the ir resonance using a self-consistent linear response theory 
which formally takes into account the mixing of the particle-hole (p-h) with the p-p channels 
in the SC state. This formalism is explained in Sec.fl[ In Sec.pTT], we present numerical results 
based on this formalism and show the overall structure of the spin correlation function. We 
then give an approximate but analytic expression for the resonance in Sec.|V|. In Sec.|V|, we 
compare our formalism with the results obtained by using equations of motion for the t-J 



model and with the SO (5) quantum non-linear sigma model. In Sec.|V]], we summarize the 
results and conclude the paper with some general remarks. Before going into these details, 
we give here some general features of the 7r-resonance. 

The central object of the theory of the 7r resonance is the so called n operator [7], defined 

by 

= ^2(cosp x - cosp y )cl +Q a a a y cl p (1) 

with a a being Pauli matrices and cj = (cL, cL) . This operator is a spin triplet and carries 
charge two. This operator inspired one of us (SCZ) to formulate the SO (5) theory of high 
T c superconductivity [p2] . Together with the total spin and total charge operators, the six 



7i operators form an 5*0(5) Lie algebra. A natural vector representation of this 50(5) Lie 
algebra is the superspin 

n = (ni,n 2 ,n 3 ,n4,n 5 ) 
= ( / ±±^,N x ,N y ,N Z} / ^^) (2) 

formed out of the antiferromagnetic (AF) order parameter 

^4^ c V% ( 3 ) 

and the real and imaginary components of the <i-wave superconducting (dSC) order param- 
eter 

A = Y,9 P C-pic P p (4) 
v 

Here, 

9 P = -^(cosp x - cosp y ) (5) 
is the d-waye form factor. The tt operator rotates N a and A into each other 

[ir a ,Nfj\=iA6 af } (6) 
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therefore within the SO (5) theory, AF and dSC are unified into a common object, called 
superspin, which can be pictured as a unit vector on an 50(5) sphere, see Fig.|I| A direct 
first-order transition between these two phases can be induced by a chemical potential [/,, 
and the superspin flops from the AF direction into the dSC direction. However, inside 
the dSC phase, there are 4 collective modes, which can be viewed as Goldstone modes 
of the spontaneous SO (5) symmetry breaking. The usual SC phase mode corresponds to 
the rotation inside the dSC plane, while there are three extra ir modes, corresponding to 
rotations towards the AF directions, see Fig.|l|. Because fj, breaks the SO (5) symmetry 
explicitly and constrains the superspin to lie at the equator, the 7r fluctuations are massive. 
From this general consideration, we expect its mass, or the resonance energy, to decrease 
with decreasing doping. 



( Neel ) 

n {2,3,4} 




FIG. 1. Geometric interpretation of the tt resonance in the superspin phase. 

The SO (5) theory therefore provides a simple picture of the it modes as collective rota- 
tions in the SO (5) sphere. Without going into the microscopic details, this picture imme- 
diately provides us with some useful quantitative information. Inside the dSC phase, the 
right hand side of the operator equation (H) can be replaced by the expectation value of the 
dSC order parameter, giving 

[tt„, N p ] = i(A) 8 aP (7) 

This equation can be simply interpreted as the commutation relation between a set of canon- 
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ically conjugate variables, just like p and q in elementary quantum mechanics. Therefore, 
we see that a new set of collective quantum degrees of freedom emerges in the broken sym- 
metry state characterized by a dSC order A. This simple consideration explains why the 7r 
resonance mode is only observed below T c . However, in the regime where a fluctuating A 
order parameter exists, the 7r resonance can also appear as a broad feature. 

In the dSC state, the n mode couples directly to the spin operator Sq = N. What is the 
dynamics associated with the 7r operator? This is a model dependent microscopic question. 
If one is dealing with an 5*0(5) symmetric microscopic model |]23|-p5l, the dynamics of the 
7r operator is determined by the equation 



5 "a 



u < (8) 



where ujq = —2fi. Therefore, 50(5) symmetric models predict a sharp 7r resonance whose 
energy scales with doping. The dynamics associated with the coupled 7r and spin operator 
in the t-J model is the central question studied in our current paper. However, even without 
detailed microscopic calculations, we can give general arguments to estimate the contribution 
from the 7r operator to the spin correlation function. Equation (|?p leads to an important 
sum rule for the mixed correlation function between the spin and the n operators. Defining 
the mixed correlation function as 

and making use of eq.(^), we have 

J ^m a p{u) = - MO|A|0) (10) 

In addition, we also have another sum rule for the ir correlation function, which follows from 
the commutation relation 



ir a , tt^ 



[l-n)5 af3 (11) 



where n is a filling factor (half filling corresponds to n = 1). From these two sum rules, we 
can put a lower bound on the tc contribution to the spin excitation spectrum as 



4 = - / ^Im X + -(Q, W ) >^L (12) 

71 J n peak 1 — Tl 

as shown in the Appendix. Therefore, one would expect the ir contribution to the dynamic 
spin correlation function to scale as the square of the dSC order parameter and inversely with 
the doping x = 1 — n. Both of these conclusions are consistent with the experimental finding 
in the optimally doped YBCO that the neutron resonance mode disappears above T c as a 
sharp excitation, and with the doping dependence of its intensity. We can use typical values 
of A( 7r ) = 40meV (see Appendix for converting the order parameter to the quasiparticle 
energy gap), J = 120meV and doping x = 15% to estimate the lower bound for I w as 0.32. 
This is close to and consistent with the experimentally measured value, 0.51 ±0.1 ||. In a 
realistic model, the it operator is not an exact eigenoperator of the Hamiltonian, and eq.(^|) 
is only approximately fulfilled with other contributions to the energy uq. However, as long as 
the 7r operator remains as an approximate eigenoperator, it will make a sharp contribution 
to the spin correlation function, and the energy of the mode will have a leading contribution 
of -2/i. 



II. SELF-CONSISTENT FORMALISM 

In this paper, we shall study the tt mode of the t-J model with nearest-neighbor (n.n.) 
hopping. Before presenting the details of the formalism, we would like to answer some 
general questions regarding the use of this model and the approximations. 



The first question concerns the effect of the n.n. Coulomb interaction V Ylij UiUj p6 
Even if we did not include a bare V term, the reduction from the on-site Hubbard model 
to the t-J model would generate such a term with V = —J/4. Actually, at this particular 
value of V, the interaction between the quasi-particles making up a spin triplet is zero. One 
might be concerned that without the multiple scattering in the triplet channel, there would 
not be any ir resonance. However, this is not the case. Even in the absence of a triplet 
interaction, there is a sharp it mode given by 
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J dte^e(t)(0\n a (t)4(0)\0) = u + 2 n +iQ <W (13) 

This occurs because the p-p continuum collapses to a point at total momentum q = Q. 
Interaction in the triplet channel simply shifts the resonance energy from — 2/i. In this 
paper, we shall only calculate the tt and spin correlation function with V = 0. The effect 
of V is two-fold, it changes both the interaction in the tt triplet channel and the energy to 
destroy a c?-wave pair in the ground state, thereby changing the chemical potential /i. Since 
the V interaction does not distinguish between the singlet Cooper pair and the triplet tt 
pair, these two contributions will essentially cancel each other. This cancelation is indeed 
observed rather accurately in the numerical calculations in both the Hubbard and the t-J 



model p0| , [2T|| . Because the J interaction is different in the singlet and triplet channels, this 
cancelation does not occur. Therefore, in this paper, we shall only study the effect of the J 
interaction. 

The second question concerns the effect of the next nearest neighbor hopping term t' 
2D| , PT|pB| , [27[] . In the presence of this term, the p-p continuum no longer collapses at total 



momentum Q, and it is not clear if the tt mode can remain sharp in the presence of t'. 
This question depends on the bandwidth around the (tt, 0) and (0, tt) points in reciprocal 
space. While the bare bandwidth might be large, it is known from both photoemission and 
numerical experiments that many-body corrections reduce the bandwidth at these points 
significantly. Assuming the reduced band structure, the tt mode remains sharp in the normal 
state. Direct numerical calculations on the tt resonance also show that the tt mode remains 
sharp for a wide range of t' pO| , pT| . Because the many-body reduction of the bandwidth 
is hard to obtain from direct perturbation theory, we shall not address the t' issue in this 
paper. 

In this work, we shall mainly discuss the two-dimensional case where the tt operator 
carries momentum (n, tt). Generalizations to bilayer system is straightforward. In this case, 
the tt operator rotates the 3D AF state into the 3D dSC state, and carries momentum 
(tt, tt, 7r), i.e. it is odd under bilayer interchange. If the 3D tt operator is an approximate 



eigenoperator of the inter-layer Hamiltonian, analysis presented in this paper will carry 
through in the bilayer case as well. 

Finally we would like to address the issue of the large Hubbard U repulsion or the no 
double occupancy constraint in the t-J model p8|j29| . In this paper, we shall only treat the 



Hubbard U within the Hartree approximation. In this case, its effect can be captured by 



a renormalization of the chemical potential |29[ and the hopping t. Alternatively, we can 



treat the t-J model within the slave boson mean field theory. Here one replaces the electron 
operator q ct by a product of bif ia . Within the dSC state, the holons bi are condensed and can 
be replaced by its c-number expectation value. The resulting Hamiltonian for the spinons 
fit, is just a t-J model with renormalized parameters, where the constraint is only treated 
on the average, again by adjusting the chemical potential and renormalizing the hopping 



parameter [pOf] . These two formalism therefore lead to the same perturbation series in the 
interaction J. 

We now review the self-consistent formalism for computing the spin correlation function 



in the SC state. This self-consistent approach has been pioneered by Anderson [31] and 



Rickayzen |52| in treating the problem of the response of a superconductor to an electro- 
magnetic field and later used by Bardasis and Schrieffer to study collective excitations in 
a superconductor |33]. The basic idea of this method is the same as that of any linear 



response calculation. We perturb a system by a small external field and then compute the 
corresponding induced response. It is however important to remember that when the system 
has SC order, any fluctuation in the p-h channel immediately mixes with fluctuations in the 
p-p, and hole-hole channels. This mixing is responsible for restoration of the transversality 
of the electromagnetic response of a superconductor J32| and preserving the Ward identities 
f36|j . Microscopically it corresponds to taking into account the response of the supercon- 
ductor due to the backflow of the condensate as well as the creation of the quasi-particle 
excitations. We have applied this formalism to the //-resonance in the negative-f/ Hubbard 
model |3B|] ( see also |37J ), and shown that it constitutes a conserving approximation, which 
gives excellent agreement with the exact theorems on the ^-resonance of the U < Hubbard 
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model p4,35|. Similar formalism has been used recently by Kohno, Normand and Fukuyama 



Salkola and Schrieffer |39| , and Brinckmann and Lee |^D| to study collective excitations. 
In this paper we emphasize that the origin of the neutron resonance peak is coupling 
to the p-p channel below T c . The SCLR formalism is a complete framework which takes 
this effect into account, and has been shown to agree with exact theorems where they are 
available [j36[| . However, the naive RPA formula Xrpa = Xbcs/{^ + VqXbcs) also contains 
partial information about mixing into the non-interacting p-p channel due to the anomalous 
F^F term in Xbcs- Therefore the peak observed at —2\x in the RPA treatment may also 
have its origin due to p-p mixing. This argument is further strengthened by the findings in 
our present work that the RPA peak at — 2fi moves to the energy of the interacting triplet 
pair within SCLR formalism. 

We start by considering the t-J model in the presence of a magnetic field Re hg^e 1 ^'^ 
( only the Zeeman effect of the applied field is of interest to us ) 

« = -tE (clc ja + h.c.) + J SiSj- - S-h^e** clc ia (14) 

<ij>a <ij> 

where S~ = J2 P c li c q+pi is the spin operator. The system responds to the applied field h qu 
in the spin channel as well as in the n channels in such a way that the operators 



Sq — ^2 c I+p1 c pI (15) 
p 
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q -E%W C -PT (I 6 ) 



J2g P c~q-pic P i (17) 



get non-vanishing time-dependent expectation values. Their Fourier transform will be de- 
noted as S quj = J dte iuJt (S q (t)) and 7T± = / dte luJt (rf (t)) . The weight function, g p = 
(cosp x — cosp y )/\^2 (as defined earlier in eq.([|)), of the operators arises from the as- 
sumed c?-wave symmetry of the order parameter^ 



In principle, the symmetry of the interaction leads to responses in the p-h channel that differ 
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The Hamiltonian ([[4]) is then linearized by factoring out the dSC gap|^ A p = A g p and 
the quantities 7r^ and S qLU : 



pa p p 

J 
1 



P 



+ £ 9pc-p,c q+p , + £r" e"^ £ ^ 



5-P-L 



(18) 



where = J(cosq x + cosq y ) and e p = — 2t(cosp ;c + cospj,) — /i. We proceed by taking the 



first two terms in eq.(18) as the unperturbed Hamiltonian Tio and use Kubo formulas to 
treat the last three terms as the perturbation TCi, 



(/(*)>=-</ dt'{ /(t),Wi(f) ) Ho - 



(19) 



This procedure is described in detail in our earlier paper on the rj excitation of the negative-C/ 
Hubbard model |JB . 



It is convenient to introduce the amplitude and phase oscillations as = 7T+ + n and 



qui qui 

and S, 



— n u . After some simple calculations, we arrive at the coupled equations for b qaJ , 



qui) 



qw 



s, 



J 



qui 



m+b 



quj 



J 

T 



m - h qu ~ VqXo(S qLU - hq u /V q ) 



(20) 
(21) 



from (|l5|) by the possible symmetry factors S q = J2 P a p c p + q -\ c pi where a p may be (sinp x ± sinp y ) 
or (cosp x ± cospy). Such 5°^ fields provide intermediate states that the spin fluctuations can be 
scattered into. In general this can modify the amplitude of the induced S qu field. However it turns 
out that for the q's of interest near Q = (tt,tt), the effect of such S quJ fields is negligible. 

2 Ao is the energy gap, which is related to the order parameter A defined earlier by the relation, 
A = Vacs A with V B cs = ^- 
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where 

dv 



t++= i Hi J ^i G -p( u - u} ) g ph{~ v ) + G- p (v)G p+q (u - v 
+ 2F_ p (v - uo)F p+q {-v)} 

= 2 Ys^iU-pUp+q + V- p V p+q ) 

dv 



\2 v pq 

/' ^ ~ V pq 



i Y,9lJ i^{G- P {v - uj)G p+q {-v) -G_ p {v)G p+q {u -v)} 
2 Y,9 2 p iy 2 _ p v 2 p+q - u 2 _ p u 2 p+q )- 



LU 2 — V 2 



dv 

% 

v 

2F. p (v - u)F p+q {-v)} 



^ g2p I 2^ G ~ P ^ ~ ^ G p+^~ u ^ + G ^ p (v)G p+q (uj - v) 



^ X/ 9p( U -p U p+q ~ V -p V p+q) 



2 "PQ 



ip\ — p-p-rq "-P"ptqj o „2 
p U ~ V pq 

dv 



m += i^dpj ^{Fp+q(v)G-p(v + u) - F-p(-u - uj)G p+q (is)} 

= 2 ^""^ 9p u P+q'Up+q( U -p V —p) o ' 2~ 
p V - V p q 

-^dpj 7^{Fp+ q {v)G- P {v + + F_ p {-v - u)G p+q {v)} 



p 



+q u p+q 2 ,,2 

p ^ - U pq 

dv 



Xo = ^ . 



E/ ^ {G- p (z/)G p+g (z/ + u) + F^ p (-v)F p+q (v - u)} 
Y,(u p+q v^ p - u. p v p+q ) 2 Ppq (22) 



In the equations above, E p = yje 2 + A 2 , i/ pg = E p+q + E_ p , u p v p = A P /2E P , u 2 p = + J^) 5 
and v 2 — ^(1 — ^-). The Green's functions have been defined as 



G P M = J dte^(-i)(Tc pa (t)cl(0)) 
F p (cu) = J dte^(-i)(T CpT (t)c_ pi (0)) 

F» = / rfte^(-0(TcL pi (t)c pT (0)) (23) 



In eg. Q2"!^) , uo should be taken to have an infinitesimal imaginary part, T = + , coming from 
causality. 
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Solution of eqs. (pTD gives the dynamical spin susceptibility in this Self-Consistent Linear 
Response Theory (SCLR). 



XSCLR(q,u) = i J 

It can be written in the form 



dt^9(t)([S}(t),SZ q (0)]) = ^-. 



XSCLR(q, v) 
Xirr 

Ax 



Xi 



1 ^qXirr 

Xo + Ax 



J mi + mi — ^m\t — ^mit ++ + ^m_m + t 



J _2 



J 2 + 2 



(24) 



(25) 
(26) 
(27) 



4"++ - 4"— -r i6"+H 

and may be understood as a modified Random Phase Approximation (RPA) where the bare 
bubble Xo has been modified by including the ladder diagrams. Fig.[2| gives the diagrammatic 
interpretation of formulas (|26|) , 





T 



+ 



T 



FIG. 2. Modification of 

Xirr due to ladder diagrams 



The procedure for finding the p-p correlation function 



P(q,t) = -i0(t)(\7r+(t),7r-(a)]) 



(28) 



is very similar to the one shown above for the spin channel. We only need to add an external 
field in the 7Tq channel and compute the response in the same channel. Skipping the laborious 
but straightforward calculations we present the final expression for its Fourier transform 



12 



P(q,u) 



t 



■+ 



+ t 




2t 



■+ 



4 



— )+m 2 _(l-|t ++ )-2m + m-(l-|t + _) 



(29) 




t++ ~ jt ~~ T6^+- + 16^++^ + J 2 



m 2 + (l-±t — )+m-(l-|t++)+f m+m-t+- 
l+V,xo 



In the normal state this reduces to a simple T-matrix expression that was studied in 
0]. There it was shown that, in the normal state, the p-p spectrum at Q is dominated 
by the collective 7r-mode resonance that appears due to the collapse of the p-p continuum 
( e p+Q + e -p — — 2/i ) and the repulsive interaction of two particles in a triplet state sitting 
on n.n. sites. We suggested that this collective mode may contribute to the spin- fluctuation 
spectrum when the system becomes superconducting. However such argument raises an 
immediate concern that superconductivity could in principle lead to another effect — a 
significant broadening of the ix resonance, due to possible scattering into the p-h excitations. 
The goal of the next part is to show that this does not happen. The ir resonance survives as 
a collective mode and affects strongly the dynamic spin-spin correlation function in the SC 
state. The important point here is that unlike xo, ^X i n Xirr contains information about the 
7r resonance. As we shall see in the next section, Im \%rr nearly vanishes at the ir resonance 
energy, where Re Xirr is sharply peaked. The combination of these two effects give rise to a 
sharp ir resonance. 



It is well known that the RPA form of the spin correlation function overestimates the 
antiferromagnetic instability. Therefore, if we see a peak in the dynamic spin correlation 
function, it is important to check if it is an artifact due to the RPA type of overestimate 
or due to a genuine collective mode. Moreover, the size of the dSC gap relative to T c is 
significantly larger than the BCS estimate. The BCS gap equation for Hq is 



with the bare pairing interaction Vbcs — 3 J/2. However, this gives a 2A( 7r! o)//c B T c ratio of 
the order of 4 which is small compared to the typically observed value of 6 to 8. Therefore, 



III. NUMERICAL RESULTS 




(30) 



13 



in what follows, we take two approaches to these problems. We introduce an effective 
reduction of the antiferromagnetic vertex Vq = aV^ re with a < 1, as a way to model vertex 
corrections, or we take the dSC gap A to be bigger than its mean field value. Both of these 
approaches have the effect of removing the RPA type of AF instability. We shall see that 
the n resonance is robust against these variations. 

A. 7r resonance and its robustness against vertex corrections 

In this section we take J = O.Qt and /x = — O.St. We choose the mean-field value of 
A = 0.0094t and the reduction of Vq is set by a = 0.82. We assume a finite value 
T = 10~H for the imaginary infinitesimal in the energy denominator and perform integration 
by dividing the Brillouin zone into a 32000 x 32000 lattice. 
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FIG. 3. Re Xirr, Im Xirr, Im XSCLR, and ImP vs. w 

This plot shows the "mechanics" of the resonance in Xsclr- in the normal state the 
p-p channel has a sharp peak at luq « — + 4(1 — n) = 0.6551 Notice that there is no 
visible shift of the energy of this resonance in the SC state, but only a small broadening. 
This resonance in the p-p channel P(Q, to) then leads to a peak in Rexirr- Consequently at 



a frequency where ReXi 



i 



the real part of the denominator in the SCLR expression 
( p5T) vanishes leading to a peak in ImxscLR- At these frequencies, the imaginary part of 
the denominator (Imxirr ) is also small, and the resonance appears to be quite sharp. 
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FIG. 4. Hexirr and Re^o vs - w - The dotted line represents the line of 1/2 J. 



In Fig.H, we compare the real part of Xirr with that of Xo- As discussed above, we have 



resonance peaks in ImxscLR when Rexir 



i 



We can see that taking x%rr instead 



of xo considerably suppresses the divergence around — 2/x (this divergence comes from the 
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dynamic nesting of the Fermi surface. It gives rise to the RPA peak, the only resonance one 
gets from a naive RPA calculation.) and leads to the development of a peak at the energy of 
the 7r-excitation. It is easily noticeable that if we do not take into account reduction of Vq, 
but exploit the bare value of Vg are = —2 J, then Rexirr will cross it at two points (u ~ — 2/x 
and lvq), giving rise to both — RPA and it peaks (see Fig.[7]). However, the divergence of 
Kexirr around ujq is much stronger, making the ir peak more robust against variations in Vq. 
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FIG. 5. Im Xirr and Im xo vs w 

In Fig.^|, we show the Imaginary part of Xirr and xo- Note that a dip develops in Imxirr 
at the energy of the 7r-excitation. This means that the 7r-resonance is much less damped than 
one might have expected. In the normal state, the stability of the ir resonance is guaranteed 
by the absence of the phase space available for decay (p-p continuum collapses to a point). 
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In the dSC state, this argument no longer works. Mixing of the p-h and p-p channels could 
provide a mechanism for the decay of the 7r excitation. However, we see that the system 
accommodates the 7r excitation by suppressing Imxi rr at its energy. In Section |IVL we shall 
give an approximate analytical derivation of this important feature. 



Im Xsclr 

lm Xrpa 




JjA I I 
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CD 

FIG. 6. Im xsclr and Im xrpa vs. lq 

In Fig.|6], we compare the self-consistent spin-spin correlation function Xsclr with the 
one obtained from the RPA calculation, Xrpa- The latter one has an RPA peak that comes 
from the dynamic nesting of the tight-binding Hamiltonian at momentum Q. In Iuixsclr, 
this peak disappears almost completely, and the the spectral weight is transfered into the 
7r-excitation. 
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FIG. 7. Im xsclr vs. to for two different values of vertex correction parameter a. 

In Fig.|7|, we show the comparison of different choices of a in Vq. Notice the coexistence 
of the RPA peak with the 7r peak for the choice of bare parameter (a = 1). Reducing a has 
no effect on the it resonance but completely destroys the RPA peak. 

From the analysis carried out in this subsection, we conclude that the RPA peak might 
be the the result of overestimating the AF instability, while the 7r peak is robust against 
vertex corrections. 



B. Robustness of the it peak against variations of the superconducting gap 

Another way of suppressing AF instability within RPA or SCLR formalism is to choose 
a larger dSC energy gap. 
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FIG. 8. Ivuxsclr vs. to for different values of Aq. 



In Fig.^, we compare the results of SCLR calculations for the spin correlation function 
for two choices of Ao- The smaller one, Ao = 0.0094i, corresponds to the self-consistent 
(mean-field) value, and the bigger one was taken as Ao = 0.05t. In these calculations, we 
take J = 0.6t and n = — 0.3t as before, but with Vq = —2 J, the bare value (a = 1). 

We observed in the previous subsection that two peaks (RPA and ir) coexist with the 
choice of the mean-field value for Ao and a bare value for Vq. Fig|| shows that taking a 
larger dSC gap removes the RPA peak and increases the spectral weight of the ir peak. This 
has an even stronger effect than we saw in the previous section by reducing the AF exchange 
constant. The latter one, as we found, only removes the RPA peak without affecting the n 
resonance. It is also interesting to find that for the larger gap there is an increase in the 
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energy of the resonance. 
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FIG. 9. Xirr and Xo vs w f° r the case when Aq is larger than its BCS value. 



A tenacious effect of the large dSC gap is explained in Fig.||. Here the choice of parame- 
ters is the same as in the previous Figure with A = 0.05t. By looking at Xirr m this case of 
large Ao, we find that the RPA peak in the real part (u « — 2/i) has completely disappeared. 
For the mean-field value of Ao, there was only a suppression of this peak. In contrast to 
that, the only effect of taking a larger A on the tt peak in Re\ irr was to make it broader. 
This broadening explains the increase in the total weight of the 7r peak in Xsclr (the slope 
of Rexirr at the crossing point with 1/2 J determines the total weight of the 7r resonance in 
Xsclr- See more on that in Section [IV}). Also note an enormous suppression of imaginary 
part of Xirr for energies below uq in this case of large A . 
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FIG. 10. Im XSCLr(<1, for q = (q, q) 

On Fig.pl], we show the q and ui dependence of Im XscLn(<i, This plot has been done 
for Ao = 0.05t, a = 1 and for computational reasons we took a larger T = 10 _2 t, which leads 
to considerable smearing. However the general picture of the q dependence of ImxscLi?(q, ^) 
may be seen quite clearly. We have an incommensurate structure at low frequencies, which 
is followed by a commensurate 7r-peak. Right after the peak, there is a missing spectral 
weight at the commensurate wavevector. Qualitatively this picture is similar to what is seen 
in inelastic neutron scattering in YBCO ||. (The incommensurate peaks in that case will 
be rotated by 45 degrees due to a different band structure. Calculations for the YBCO band 
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structure will be published elsewhere ffl| . ) 



IV. ANALYTICAL DERIVATION OF THE vr RESONANCE 

In the previous section, we showed numerically that at the energy of the 7r-excitation, 
Xirr possesses a sharp peak in the real part and a dip in the imaginary part. In this section, 
we study the origin of these properties analytically 3 . When looking at the structure of the 
expressions in eq . (|2lD , one encounters very often analogous integrals that only differ by a 
factor in the p-summation. To simplify the following analysis, we make an approximation 
that this factor may be replaced by its average value of 1. A similar assumption has been 



used in [22] to obtain approximate SO (5) algebra. It is important to realize that one should 
take the average of gi not over the whole Brillouin zone but over a narrow band around the 
Fermi surface, since in most of these expressions the other factors in the integrals restrict 
the important domain of integration to this region. 
We introduce 

t t. ,\ _ v r? 1 ~ v -p~ v p+Q 
P V ~ v vQ 

J( \ \ - n Up+QVp+Q - U-pV^ p 

Wu)-2^9p — 2 — "2 K 6V ) 

v u v vQ 

and use some identities for the BCS coherence factors and the approximation ~ 1 to 
express all the factors in (^) as 

t++ = -Afil^u) 

t + _ = — 2uli(u) 

1 — n u 2 T . . 
= JiH 

m + = —2fil 2 (uj) 

m_ = —ujI 2 {oj) 



3 We restrict this analysis to the case q = Q = (it, 7r). 
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2 uj 2 - 4u 2 T . . 



Substituting these expressions into eq. (|27|) and using the identity — 4/iA /2(^) = (1 — 
n) + (4/i 2 — uj 2 )Ii{uj) which holds within the approximation described above, we get 



v 1 , jgM («>* - uftju,* - 4f) 

V " T " Vbcs 2A uj 2 -uj 2 + A JI 2 (uj)(uj 2 -2^ ) 



where Vbcs comes from the gap equation J2 P 9p u p v p — ^o/Vbcs (eq.(|30|)), and co> = —2// + 
^(1 — n). In the mean field analysis of the t-J model, Vbcs — 3 J/2. Since < for 

uj ~ ujq, the denominator of the expression (|33"D vanishes when the frequency is larger than 
co>o but very close to it (the energy separation is proportional to Aq ). This explains the peak 
in KeXirr, and the factor uj 2 — uj 2 . in the numerator explains the dip in the imaginary part. 

Expression ( |33| ) allows us to estimate the integrated spectral weight of the 7r-excitation 
in the spin-spin correlation function. For simplicity, we neglect a small imaginary part of 
Xirr near ujq. Then a pole in Xsclr occurs when 

ReXirr = ~^r- (34) 
Expanding xsclr around this frequency cu* we find that 

Xirr 



Xsclr 



V c 



! = ' =*_. (35) 



Earlier we introduced two distinct energies in our system. The BCS coupling Vbcs and 
the AF coupling Vq = —2.1. If we consider a hypothetical situation when Vq = —Vbcs, 
we can see that the condition (K3j) is satisfied exactly at uj and a simple calculation gives 
Xsclr(oo) = v ^l n) u-Zo+io - If we take nere V Q to be V bcs = 3 J/2, we have for the intensity 
of the 7r resonance as 

! '=nL *^=wfrj- (36) 
The right hand side is equal to the expression derived in the Appendix as a lower 
bound. In eq . (p6|) , v and v' characterize the width of the 7r resonance around ujq, 
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and we introduced a factor 1/tt since the Lehmann representation of ([24] ) is given by 

2 

Im x(Q,u) = ^J2n ( n \So\fy S(co — u} n0 ). This definition of I n is the same as A(T) of 
ref. Expression (|36| ) is also what we obtained for the intensity of the 7r resonance in |7j 
using a T-matrix analysis^. 

If we take realistic values of J = 120 meV, A( 7r ) — 40 meV (this corresponds to 
A = 28.3 meV) and 1 — n — 15%, and substitute them into eq. (|3~6D , we get I n = 0.32. 
For the t-J model, \Vq\ > Vbcs, and the energy satisfying the condition ( |3~4"|) is lower than 
ujq and the slope 9 ^ r is smaller than the above estimate (see Fig.|j). This will be partly 
cancelled with the increse of Vq, and we expect eq. (|36|) gives a semi-quantitative estimate 
for 

V. COMPARISON WITH THE SO (5) EQUATIONS OF MOTION 

We now study the Heisenberg equations of motion (EOM) for the 7r and spin operators 

n pQ = c p+QT c -pT 

^pQ = C -P-QI C PI 

SpQ — C p+Q] C pl (37) 



using ([Til) as a Hamiltonian. A closed set of equations may be obtained by taking com- 



mutators of the operators (|37|) with the Hamiltonian and then factorizing the results in 
terms of the occupation numbers for the electrons = (c^Cpa) and BCS anomalous aver- 



ages UpVp = (c_pfC p |). As shown by Anderson and others f31| , this procedure recovers the 
Modified Random Phase and T-Matrix Approximations. 



4 For unit convention, see Appendix B of ref. 0. 

5 Note a factor of ^/2 difference in the definitions of Aq here and in ref. Aq (present) 



v / 2A (ref. 
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*Z) + i (! - 



-p) n pQ 



3J 



+ S+p-g q) I] ujwCp - A;) + 4 J Mp^p 51 S, 



kQ 



(38) 



n,% 



vQ 



kQ 

k 

3J J 

z fc z fc 

+ e-j>pQ - g ( X - U J - Z n k Q v(p ~ k ) 



3J 



+ -7^( S pQ ~ S - P -Q q) Z u k v k r](p -k)-AJ u p v p S, 



kQ 



(40) 

* p fc 

Here ?y(p) = cosp x +cosp y and the bare dispersion is renormalized into e p = e p — ^ Yk v \ v(p~ 
k). The latter corresponds to a trivial rescaling of t which we will disregard. 

The operator of the collective excitation in the ir + channel is ttq = YpQip) 71 ^- Then 
from eq.(^), we have 

JA 



^Q + 77 b Q 

VBCS 



(41) 



where Sq = J2 P SpQ an d Vbcs has been defined earlier. Analogously, 

JA 



H, 7r; 



-^OTTQ - - 



-Or 



BCS 



(42) 



We can see that in the SC state the EOM for ttq no longer close on themselves. The 
third and the forth terms in eqs.([38"l) and flip]) , that come from anomalous self-energy and 
scattering correspondingly, do not cancel each other exactly. This may be contrasted to the 



77-excitation in the negative-f/ Hubbard model [36], where exact cancelation of such terms 
occurs. 

The collective mode in the S channel may be obtained by summing eq.([3^) over different 

p's. 

2JA 



V, 



BCS 



7Tn — 7IV 



(43) 



In order to derive this result, we had to disregarded the first term of eq.(^). In the language 
of our earlier SCLR approach, this means neglecting xo i n comparison with A%. In a close 
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vicinity of uj , this is a justifiable assumption, because at these frequencies Ax is strongly- 
peaked and is the dominant part of Xirr- However, it is less so at other frequencies, where 
the incoherent continuum is more important. Thus the meaning of going from eq.(|39"D to 
eq.(^3|) is a single mode approximation, which captures collective degrees of freedom only. 
How good is this approximation? One can see from the numerical results of Sec.^TJ that for 
Ao around 0.1 J the 7r-peak already became a dominant feature of the Sq spectrum. Some 
estimates of the realistic value of Ao find it to be close to 0.2 J. In this case, such single 
mode approximation will truly be a good one []. 

It is instructive to compare our microscopic EOM's with the SO (5) EOM's in the SC 



state 22 . We have 



-in r 



±(B l5 n± + g(n 5 )n a ) (45) 
--L(n 5 )U-7r a ) (46) 

with = Li a ±iL 5a and we assumed dSC ordering along n 5 . Equation (|4"5| ) is the analogue 
of eqs.(|4]) and (|42"D, and eq. (f4~6"|) corresponds to eq. (^3|) . Results in this section therefore 
give a microscopic justification and quantitative derivation of the 5*0(5) quantum nonlinear 
cr-model as a long- wave- length theory in the dSC state [E2j. 



6 However, this approximation seems to overemphasize the importance of the anomalous scattering 
terms in the energy of the n excitation. The anomalous self-energy tends to increase the energy 
of the resonance, whereas the anomalous scattering decreases it (see eqs.(FS^), p9|), and (fiC|), for 



example). In Sees. Ill and IV, we saw that in the complete calculations, the resonance energy in 
the SC state turns out to be above its value in the normal state. However, from eqs.(|4l|), ( f42| ) and 
|), we find the resonance energy in the SC state to be decreased 

/2JA \ 2 



2 2 
UJ S = u 



[v^) ■ (44) 
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VI. SUMMARY 



We have presented detailed analytical and numerical calculations for the contribution 
from the ir resonance to the spin correlation function in the dSC state. The results of these 
calculations support our earlier interpretation of the resonant neutron scattering peak in 
terms of the p-p collective mode in the 7r channel. Various approximations were used in 
the calculations presented in this work, some of them are model-dependent and may not be 
well-controlled. Therefore, it is important to summarize here the main points leading to our 
conclusion. 

1) From general model- independent sum rules on various correlation functions, one can 
conclude that the contribution from the ir correlation function to the dynamic spin cor- 
relation function is of the order of | A| 2 /(l — n), in excellent agreement with the two key 
experimental observations, namely the vanishing of the sharp mode above T c and the doping 
dependence of its intensity. 

2) Within model-dependent calculations, there is a well-defined ir mode in the p-p channel 
in the normal state, and this mode couples to the p-h spin channel in the dSC state, where 
it remains as a sharp excitation. The energy of this mode is not directly related to the dSC 
gap, but is directly related to the doping x. In the underdoped materials, the dSC gap 
increases slightly as doping is reduced, while the neutron resonance peak energy decreases 
with x. This important experimental finding shows that the neutron resonance peak is not 
simply a "2A" phenomenon, and our interpretation in terms of the 7r resonance naturally 
resolves this apparent paradox. The doping dependence of both energy and intensity of 
the neutron resonance peak were predicted |7|,|22[| before the experiments in the underdoped 
superconductors were carried out |H|. 

3) Many approximations within our current calculations are not completely controlled. 
However, the main behavior can be verified in the case where exact knowledge is available. 
First of all, detailed exact diagonalization studies have been carried out both for the t-J 



and the Hubbard models pT|,pO|. It is clearly seen that the it mode in the p-p channel 
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exists in all doping range, and it has a low-energy peak where both the energy and intensity 
scale with x, in agreement with ourT-matrix calculation. In contrast to the 7r correlation 
function, the spin correlation function does not have sharp peaks in the high doping range. 
In the doping range where there are dSC fluctuations, the 7r peak coincides with the spin 
peak. From these results, one can conclude that the 7r mode is a genuine collective mode. 
We can also compare our approximations with the exact £0(5) models [ft3fl , where the ir 
operators are exact eigenoperators of the Hamiltonian. The manipulations presented in this 
work lead to results consistent with the exact SO (5) Ward identities. The 7r resonance is 
an exact excitation of the 5*0(5) models, and it has exactly the same doping dependence of 
the mode energy and intensity as obtained here. 

4) The distinction between the "RPA peak" and the u ir peak" in the dSC state will 
be a model-dependent one. In the dSC state, they share the same quantum numbers, and 
both are based on approximate calculations. The origin of the RPA peak may be related to 
the overestimate of the magnetic instability and we see that it may not be robust against 
variations of the vertex corrections or variations of the gap, both of which diminish AF 
instability. On the other hand, the SCLR treatment of the 7r peak is more robust against 
these variations. One can test these two approximate schemes within the exact 5*0(5) 
models. Only SCLR treatment including the 'V process agree with the exact answer in 
this case. Therefore, calculations including the 'V process is a better approximation than 
the simple RPA calculation. 

5) Within our approximations, the spin spectrum consists of an incommensurate struc- 
ture at low frequencies, a sharp commensurate peak arising from the triplet excitation in 
the p-p channel (the ir peak), and a missing spectral weight at commensurate wavevector at 
higher frequencies. These features are in overall agreement with experiments. The predicted 
weight of the it resonance agrees quantitatively with experiments. 

Therefore, while each of the above arguments are not complete on their own, the combi- 
nation of them makes a strong overall case. From the interpretation of the neutron resonance 
peak in terms of the 7r mode, we hope to learn a general principle, rather than a specific 
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model for fitting a specific experiment. In strongly correlated systems, most degrees of 
freedoms are strongly coupled, and most spectra are incoherent. Usually, only a symmetry 
principle can forbid the decay of a collective excitation. In the case of the resonant neutron 
scattering peak, we believe that it is the SO (5) symmetry principle at work, and the n mode 
is the pseudo Goldstone boson associated with this spontaneous symmetry breaking. In this 
paper, we have shown that such an interpretation is consistent with the key experimental 
facts, but it may not be the only possible interpretation. Its utility lies in the simplicity and 
generality of the principle, which can be applied to other related experiments and lead to 
new experimental predictions. 
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APPENDIX A: MODEL-INDEPENDENT ESTIMATE OF THE vr 
CONTRIBUTION TO THE SPIN SUSCEPTIBILITY 

A spectral function for any two operators A and B is defined as 



where D r / B (u) and Df v B (uj) are retarded and advanced response functions, respectively, 
|n)'s are eigenstates of the system Hamiltonian with energy E n , and |0) is the ground state 
with energy Eq. If we restrict the above summation only to intermediate states that have 
nonzero overlap with 7rt |0) , then such a quantity 



= J2[(0\A\n}(n\B\0)5(uj + E -E n ) 



n 



-(0\B\n)(n\A\0)5(u; - E + E n )], 



(Al) 




Y, (0\A\n)(n\B\0)5(u; + E -E n ) 



n:{0\n a \n)^0 
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-{0\B\n){n\A\0)5(uj- E + E n )\ (A2) 

may be regarded as the contribution of the n excitation to the full spectrum pa,b(w). We 
can introduce cu-integrated spectral weight as 



W2 
Ull 



IpXbZI = / dupXs^)- (A3) 



All of the above spectral functions (or spectral weight) are bilinear with respect to A and 
B, and have a property, Pa,a^ > for u > or uj 2 > oj\ > 0. Therefore, the Cauchy-Schwarz 
inequality holds provided that the same frequency condition is satisfied] 7 ] : 

PA,A^PB,Bt > \pA.m? ( A4 ) 



Here p A ,Bt can be either of p A ,Bi(u), P^ B t( w ) or IPabi. 



U12 



In Section | we derived two sum rules / duip* 01 t (u) = 1 — n and / duip* a , rt (^) = «A. If 

n a ,n a ir a ,r\ a 

most of the 7r spectrum is accommodated in an interval (uq — u, ujq + u') on the positive real 
axis and around the 7r-resonance energy ujq, we can write 



P na t 



UJQ+V 1 



~ 1 — n, 

u>o—v 



p na 



~ iA. (A5) 



Equation (|A4|) then immediately gives us 



7 This may be proved as follows. Prom bilinearity and (semi-) positivity, we have 

< PA+XB,At+X*Bi 



= Pa,ai + Vs.At + **Pa,bi + |A| p fliB t 

for any complex number A. Defining 9 as the phase of the 'mixed correlation function' as p A B f 
p* B A x = \pa,b\ 1 an d choosing A = xe ld (— oo < x < oo), we have 

p AAt + 2x\p ABt \ + x 2 p B m > 

Since this inequality holds for any real number x, the inequality ([A3) should hold. 
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p 



u) +i/ |A| 2 
>- — 



n 



(A6) 



The left hand side of this equation represents the contribution of the 7r mode to the spin 
excitation spectrum (e.g., \m.x zz /'x), and the right hand side gives its lower bound 0. This 
is a model-independent result. 

Noting that Inr\; +_ = 21mx zz , we obtain eq. ( 12 ) . When applying this result to the 
present analysis of the t-J model, where A = Vbcs^ — ^fA, we have /„- > gjs--^ 2 ^-. 

Analysis in this Appendix can be generalized to finite temperatures by considering the 
spectral function, 



Pa,b{u) 



-PE n 



- e~ t3Em )(n\A\m)(m\B\n)5(uj + E n - E, n 



(A7) 
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